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We compute the contribution of phonons to the shear viscosity rj in superfluid neutron stars, 
assuming neutron pairing in a ^So channel. We use a Boltzmann equation amended by a collision 
term that takes into account the binary collisions of phonons. We use effective field theory techniques 
to extract the phonon scattering rates, written as a function of the equation of state (EoS) of the 
system. Our formulation is rather general, and can be used to extract the shear viscosity due to 
binary collisions of phonons for other superfluids, such as the cold Fermi gas in the unitarity limit. 
' We find that rj oc l/T"", the proportionality factor depending on the EoS of the system. Our results 

, indicate that the phonon contribution to 77 cannot be ignored and might have relevant effects in the 

dynamics of the different oscillation modes of the star. 
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I. INTRODUCTION 

X/^ i After Migdal's observation that superfluidity of neutron matter may occur in the core of compact stars this 
(-H ' possibility has been seriously explored and taken into account to explain very different neutron star phenomena, such 
O ^, as pulsar glitches or the star cooling. Superfluidity would also affect many other aspects of the neutron star dynamics, 
the corresponding hydrodynamics being essentially different from that of a normal fluid In particular, rotational 

and vibrational properties of the star, or the dynamics of the star's oscillations would be specially sensitive to the 
presence of a superfluid. 

^ \ Neutron matter superfluidity occurs at low temperatures after the appearance of a quantum condensate, associated 
to neutron pairing. The condensate spontaneously breaks the global ?7(1) symmetry associated to baryon number 
conservation. In such a case Goldstone's theorem predicts the existence of a low energy mode that at sufficiently low 
^ , momentum has a linear dispersion law, and which is essential to explain the property of superfluidity. We will refer 
' generically to this mode as superfluid phonon, or phonon for simplicity. Although we are here mainly interested in 
■ neutron matter, the appearance of a phonon is however a generic feature that happens in all superfluid systems. 
• The phonons dominate the physics in the superfluid at long wavelength. At very low temperatures they also demi- 
es I nate the thermal corrections to the thermodynamical and hydrodynamical properties of the superfluid. Surprisingly, 
. even if this fact is very well-known, the presence and effects of the phonons have been largely ignored in the compu- 
tation of all the transport coefficients that are needed to study the dynamics of superfluid neutron star's oscillations 
We are only aware of Ref. 0, which studies the thermal conductivity due to the superfluid phonons in the 
inner crust of the star. 

It is our intention to compute the phonon contribution to the different transport coefficients of superffuid neutron 
matter. In this manuscript we start with the computation of the phonon contribution to the shear viscosity rj in the 
. , core of superffuid neutron stars, assuming that the neutrons pair in a ^S'o channel. We leave the case of neutron 
' pairing in other channels, the computation valid for the inner crust of the star, or the computation of other transport 
coefficients for future studies. 

In order to compute the phonon contribution to the different transport coefficients one needs to assess the relevant 
phonon collisions which are responsible for the transport phenomena under discussion. There has been a lot of progress 
in this direction in recent years. The leading phonon interactions can be completely determined by the equation of 
state (EoS) of the system, as can be seen after writing the corresponding effective field theory at leading order in a 
derivative expansion [^, Q . 

Our computation of the shear viscosity is very similar to the same sort of computations for other superffuid systems 
[l0l - [T2j . In this paper we however exploit the universal character of the effective field theory at leading order to present 
a very general formulation, that only depends on the EoS of the system. This will allow us to obtain immediately the 
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viscosity due to the binary collisions of phonons of other superfluid systems, such as, for example, the cold Fermi atoms 
in the unitarity limit [Til [l3|. Our formulation also will allow us to improve the value of 77, if a better determination 
of the EoS for the nuclear matter present in the star is achieved. 

This paper is structured as follows. In Sec. |TT]we review the effective field theory that describes the phonon self- 
interactions in a superfluid, and write all the phonon self-couplings in a model independent way, that is, without 
specifying the system we are considering. We use this theory to extract the scattering amplitudes of the binary 
collisions of phonons in Subsec. IIIBI In Sec. IIIII we present and review the specific equation of state for /3-stable 
nuclear matter that will use in our computation of the shear viscosity. In Subsec. IIV Al we use the transport equation 
to show how the value of rj can be extracted for a generic superfiuid. In Subsec. flVBI we present numerical results for 
the value of 77 when we consider the EoS corresponding to /3-stable nuclear, and also that of a free neutron gas. In 
Subsect. IIV CI we compare our results with those obtained for ^He and for the cold Fermi gas in the unitarity limit. 
We end up in Sec. |V]with a discussion of our results. In all the manuscript we use natural units, except in the plots 
when we give numerical values of different physical quantities. 



II. SUPERFLUID PHONON INTERACTIONS 



In this Section we first review the effective field theory (EFT) that at leading order in a derivative expansion 
describe the phonon interactions. We parametrize all phonon self-couplings in terms of the density, speed of sound, 
and derivatives of the speed of sound with respect to the density in a model independent way. In Subsect. Ill Bl we 
use this very general formulation to extract the scattering amplitudes associated to the binary collisions of phonons. 



A. The effective field theory 



The superfluid phonon is the Goldstone mode associated to the spontaneous symmetry breaking of a U{1) symmetry, 
which corresponds to particle number conservation. EFT techniques can be used to write down the effective Lagrangian 
associated to the superfiuid phonon. The effective Lagrangian is then presented as an expansion in derivatives of the 
Goldstone field, the terms of this expansion being restricted by symmetry considerations. The coefficients of the 
Lagrangian can be in principle computed from the microscopic theory, through a standard matching procedure, and 
thus they depend on the short range physics of the system under consideration. 

It has been known for a while that the Hamiltonian of the Goldstone mode of a superfluid system is entirely fixed 
by the equation of state [13, [l^ . In some more recent publications [1, [^ it has been realized that at lowest order in 
a derivative expansion the Lagrangian reads [^ 

Cw = PiX) , (1) 
where P{fJ.) and fi are the pressure and chemical potential, respectively, of the superfluid at T = 0, and 

X^.-a.^-^^, (2) 

where ip is the phonon field, and m is the mass of the particles that condense. After a Legendre transform, one can 
get the associated Hamiltonian, which has the same form as the one used by Landau to obtain the self-interactions 
of the phonons of ^He [3, [^ . 

This effective field theory has been used for the computation of different phonon effects in different superfluid 
systems. In particular, it has been shown [l6j that it allows to reproduce the leading thermal corrections to the speed 
of sound at low temperature that were obtained by Andreev and Khalatnikov for superfluid ^He [l7j . 

It is convenient to re-express the phonon Lagrangian in a different way. After a Taylor expansion of the pressure, 
and rescaling of the phonon fleld to have a canonically normalized kinetic term, one can write the Lagrangian for the 
phonon fleld as 



(3) 



The different phonon self-couplings of Eq. ([3]) can be expressed as different ratios of derivatives of the pressure with 
respect to the chemical potential [l6| . But for comparisons with the condensed matter literature on superfluidity, and 
for the purposes of computation in this paper, it turns out to be more convenient to express them in terms of the 
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density, the speed of sound at T = 0, and derivatives of the speed of sound with respect to the density. First of ah, 
it is easy to check that the phonon velocity is the speed of sound at T = 




'ph = J^ = Cs , (4) 



where p is the mass density, related to the particle density p as p ~ mp. 

One then finds that the dispersion law obtained from this Lagrangian at tree level is exactly Ep = c^p. The three 
phonon self-coupling constants can be expressed as (see Appendix A of Ref. 

1-2-— , 'ng= r— , (5) 



6y/mp Cs \ Cs dp J ' 6y/mp g 

while the four phonon coupling constants are 




(6) 

It is also possible to construct the Lagrangian for the superfluid phonon field at next to leading order (NLO). 
However, for our present purposes, we will not need it. It is only convenient to remember that at NLO the phonon 
dispersion law receives corrections that take the form 

Ep ^Csp + Bp^ + ■■■ (7) 

where the parameter B can be computed by a matching procedure with the microscopic theory. The value of B is 
relevant in order to see whether the decay of a phonon into two is kinematically allowed or not. Only dispersion laws 
that curve upward can allow such processes. The possibility of having these phonon decay processes is important, as 
they affect the value of the different transport coefficients of the superfluid. For the system of our interest, superfluid 
neutron matter, we have checked that B < 0, and thus we will only consider binary collisions of phonons in the present 
manuscript. 



B. The scattering cross sections of binary collisions of phonons 



A necessary ingredient for the computation of the shear viscosity is the knowledge of the scattering cross section 
associated to the collisions responsible for that dissipation. It is possible to extract them, at leading order in a 
momentum expansion, from the phonon Lagrangian Eq. ([3|). While these scattering amplitudes have been written 
down for some speciflc superfiuids (that is, specifying the value of the EoS and thus of the phonon couplings), we 
write them down here in a much more general form, that is, without specifying the EoS of the system. 
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The four Feynman diagrams of the binary collisions under consideration are depicted in Fig. [TJ The contact 
interaction amplitude is 

iMc.t. = -i\[24EpEkEp,Ek' - {EpEk p' • k' + EpEp, k • k' + EpEk' p' • k + Ep,Ek p • k' + Ep,Ek' p • k 

+ EkEk'P ■ p') + 8VX.2 (p • k p' • k' + p • p' k • k' + p • k' p' • k) } , (8) 
while the s-channel amplitude reads 

tMs = ~WG{P + K) {EpK' + EkP^- + 2{Ep + Ek)P ■ K] {Ep'K''' + Ek'P'^ + 2{Ek' + Ep')P' ■ K'] , (9) 

where we used the notation P ■ K = EpEk — 3?7gP • k, and G is the phonon propagator. The t- and w-channel 
amplitudes can be obtained from the s-channel one by using the crossing symmetry iMt = iMs{K -H- —P') and 
iMu=iMs{K ^ -K'). 

As found out in Ref. [loj . if one uses the bare phonon propagator as obtained from the Lagrangian Eq. ([3]), one 
finds collinear singularities in small angle collisions where one virtual phonon is exchanged. In order to avoid coUinear 
singularities in the intermediate steps of the computation of the shear viscosity one has to regulate the phonon 
propagator. One possibility is to resum the phonon propagator including the one-loop thermal damping. Thus, one 
writes 

^'^^^ " S2-c2p2 + iinin(P) ■ ^^^^ 

The imaginary part of the self-energy is extracted after a one-loop computation (see Ref. [iBl)- For space-like momenta 
andp<T 



imn(P) = f 1 + ^e(cy - pD , (11) 
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15c5 p \" ' Cs dp j CsP 
while for time-like momenta, and po > 0,p — > one finds [ll| 

■">n(P).Jfl(i + |)"pg5^, (12) 

Another possible way of regulating the phonon propagator is to work with the EFT at next to leading order, using 
the non- linear dispersion law Eq. ([T])- However, for the computation of the shear viscosity in superfiuid neutron matter 
wc have checked that any of these two possible regulations give the same result. This is due to the fact that the shear 
viscosity is dominated by the large angle collisions, the viscosity being insensitive to collinear collisions. 



III. EQUATION OF STATE OF SUPERFLUID NEUTRON MATTER 

In order to obtain the speed of sound at T = as well as the different phonon self-couplings one needs to determine 
the EoS for neutron matter in neutron stars. Then one can get the speed of sound, and all the derivatives of the speed 
of sound with respect to the density. A common benchmark for a nucleonic equation of state is the one obtained 
by Akmal, Pandharipande and Ravenhall [l^ (APR for short). These authors studied neutron matter as well as 
symmetric and /3-stable nuclear matter. Later on Heiselberg and Hjorth- Jensen [23| parametrized the APR EoS of 
nuclear matter in a simple form, which will subsequently be used in this manuscript. It should be remarked that the 
effect of neutron pairing is not considered in the previous references. However, because neutron pairing only affects 
those neutrons which are closed to their Fermi surface, one does not expect that the effect of pairing might have a big 
impact in the EoS. For this reason, we will neglect the effect, and consider that it only represents a small correction. 

The binding energy per nucleon (E/A) in nuclear matter consists of a compressional term and a symmetry term: 

E/A = Ecompip) + S{p)il - 2xpf = + Sou''il - 2xpf. (13) 

Here u is the ratio of the nucleon density (p) to nuclear saturation density (po = 0.16 fm^"^), u = p/po, and Xp = Pp/po 
is the proton fraction. The nucleon density is given by 



(14) 
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with pf being the Fermi momentum and v the degeneracy factor. In nuclear matter v is 4. The compressional term 
is parametrized by a simple form which reproduces the saturation density, binding energy and compressibility. The 
binding energy per nucleon at saturation density excluding Coulomb energies is Eq = 15.8 MeV and the parameter 
S = 0.2 was determined by fitting the energy per nucleon at high density to the EoS of Akmal [l^ with three-body 
forces and boost corrections, but taking the corrected values from Table 6 of [l9j . The corresponding compressibility 
is Kq = 18 £o/(1 + (5) ~ 200 MeV, within the experimental values. For the symmetry term Heiselberg et al. obtained 
So = 32 MeV and 7 = 0.6 for the best fit. 
The EoS is then 

£{p,Xp) = im + E/A{p,Xp)) , (15) 

with m being the mass of the nucleon, while the corresponding nucleonic energy density is 

eN{p,Xp) = £{p,Xp)p . (16) 

For neutron star matter made of neutrons, protons and electrons, the total energy density is the sum of the nucleonic 
contribution (neutrons and protons), e^v, and the one for electrons, ee, 

e{p,Xp,Pe) = eN{p,Xp) + ee{pe). (17) 

The pressure includes also both contributions 

P{p,Xp,Pe) =^ Pn{p,Xp) + Peipe) , (18) 

where the nucleonic and electronic contributions to the pressure are given by 

Pn{p,Xp) p,n{p,Xp) {1 - Xp) p + Pp{p,Xp) XpP - eNip,Xp), 

Pe{pe) = Me(Pe)pe - ee(Pe) , (19) 

being pi the chemical potential of each specie. Those chemical potentials are calculated as 

deN{p,Xp) _ de^jp, Xp) deN{p,Xp) Xp 



^^n{P,Xp) 
fJ-p{Pi Xp 



^ dpn dp dxp p 

dcNjp, Xp) _ dej^jp, Xp) deN{p,Xp) (1 - Xp) 



dpp dp dxp p 

Me(Pe) = yplTT^-pf, - (3^Ve)'/'. (20) 

Moreover, matter in neutron stars is thought to be in /3-equilibrium against weak decay processes among the different 
species. The constraints imposed by chemical equilibrium and charge neutrality for matter made of neutrons, protons 
and electrons are 

P"n pp ^" P'e ^ 

Pp = Pe ■ (21) 

Using the definitions for the chemical potentials of Eq. pO|) and the /3-equilibrium condition, one finds that 

Me = Mn-Mp = --^^^^^^^=4W(l-2.Tp) . (22) 

Then, one can use the charge neutrality condition for ultrarrclativistic electrons: 

p, = Ep^ = pp^ = (3^2^e)'^' = (37r2pj,)i/3 . (23) 

Using that pp ~ Xpp, wc get 

. (4^o--(^-2^.))^ . (24) 

2>TT'^ P 

From this relation, we can find an analytic equation for the proton fraction. Defining 

a = 2(45'ou^)V(^V) , (25) 
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FIG. 2: Ratio of the speed of sound with respect to the speed of light c obtained for a free gas of neutrons and for /3-stable 
nuclear matter as a function of the density. 



Eq. ([M]) reduces to 



which sets the relation between the proton fraction and the total nucleonic density. 

In our calculations for the shear viscosity we will use the APR EoS of /3-stable nuclear matter, in the paramctrization 
just presented. We have also compared our results with those using the APR EoS of interacting neutron matter and 
the EoS of a free neutron gas, that is, assuming that the interactions are very weak and their effect is neglected in 
the EoS. We have checked that the scenario of only interacting neutrons leads to very similar results for the speed 
of sound and all phonon self-couplings than that of /3-stablc nuclear matter. Those results were, however, expected 
due to the low proportion of protons in /3-stable matter. On the other hand, the EoS of the free neutron gas shows 
a different behavior to /3-stable matter, and thus of the speed of sound and phonon self-couplings, and ultimately, of 
the shear viscosity, as we will show in the following Section. 

In Fig. [2] we plot the ratio of the speed of sound Cs with respect to the speed of light c for both a free gas of 
neutrons and for /3-stable nuclear matter as a function of the density. We observe that the assumption that the APR 
model for /3-stable nuclear matter is non-relativistic breaks down at densities of the order of 1.5-2 po. For those 
densities, relativistic effects appear as the APR EoS includes not only two nucleon but also three nucleon interactions. 
This seems to suggest that at very high densities we should employ the relativistic version of the phonon effective 
field theory Q . This would imply very minor changes in the computation of 77 (see Ref. jTo| for the treatment of a 
relativistic superfluid) . 

In the following, we will perform the following approximations for the computation of the speed of sound and all 
three and four phonon self-couplings. First, to compute the mass density we only take into account the nucleonic 
part, as m 3> m,,. Further, in /3-equilibrated matter, p sa /9„. Thus, we will assume that the speed of sound can be 
computed as 



and, for the phonon self-couplings, we will substitute p in Eas. (j5l6|) . 

IV. SHEAR VISCOSITY DUE TO THE BINARY COLLISIONS OF SUPERFLUID PHONONS 



ax^ 4- 3a; — 3 = , 
with X = 1 — 2xp. This equation has one real solution given by 



(26) 



X = (-2 2i/3a + (6a2 + 2^a^{'i + 9a)f^^)/{2a (Sa' -f- ^a^{4 + 9a)y/^), 



(27) 




(28) 



We report here all the computations required to extract the shear viscosity due to the binary collisions of phonons 
in a non-relativistic superfluid. Our final goal is to extract the value of the viscosity in /3-stablc neutron matter. 
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using the EoS just presented in the previous Section. Because to leading order in a momentum expansion the main 
interactions leading to the value of the viscosity are set by the EoS of the system, we first present in Subsect. IIV Al 
a very general formulation of the computation. We follow here the same methodology used for the study of some 
specific superfluids [lol - [l^ . but we are aiming here to have a broad and general formulation that can be used for 
different superfiuid systems, or for future improvements of the knowledge of the EoS of the superfluid neutron matter 
contained in the star. 

In Subsect. IIV Bl we present the numerical results of the shear viscosity when the APS EoS discussed in the previous 
Section is chosen. To see the effect of the nucleon interactions, we also compute the viscosity with the EoS of a free 
neutron gas. In Subsect. lIVCI wc compare our results with those obtained for ^Hc and for the cold Fermi gas in the 
unitarity limit. 



A. General formulation 



The shear viscosity rj enters as a dissipative term in the energy-momentum tensor Tij. For small deviations from 
equilibrium one finds 

= -TjV,, = -77 (^d,Vj + d,V, - ^%V-V^ , (29) 

where V is the three fiuid velocity of the normal component of the system. 

The superfiuid phonon contribution to the energy-momentum of the system can be expressed as 

where / is the phonon distribution function. The distribution function obeys the Boltzmann equation Q 

where we have assumed to be in the superfiuid rest frame, and C[f] is the collision term. For the computation of the 
shear viscosity it is enough to consider binary collisions, and thus 

cm = — f , , , (2^)^^(^)(P + if - P' - K')\m\^B , (32) 

^•'^ 2Ep J {27T)^2Ek{2n)^2Ep, {2nf2Ek'^ ' ^ ^2' ' ' ^ ' 

where 

D = /(p')/(fc')(l + + fm - f{p)f{k){l + /(p'))(l + fik')) , (33) 

and Ai is the 2 o 2 scattering matrix 

M^Mc.t.+Ms + Mt + Mu ■ (34) 

The explicit expressions for the contact interaction, and s-, t- and u- channel amplitudes were given in Subsec. IlIBI 
This collision term has the property that it vanishes when evaluated using the phonon equilibrium distribution 
function 

so that C[/oq] = 0. 

For the computation of the transport coefficients we have to consider small departures from equilibrium so that the 
distribution function can be written as 

/ = ./cq + Sf , (36) 

and linearize the transport equation in 6f. For the computations of the shear viscosity, one assumes a particular 
dependence of Sf. More particularly 

6f = -hip)p,,Vki^^^^^^^^^ , (37) 
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where h{p) is an unknown function, Vki is the tensor defined in Eq. (|29|) . and we have defined the tensor 

Pki = PkPi - ^SkiP^ ■ (38) 

One can extract the value of the shear viscosity computing the contribution that this distribution function has in 
the energy- momentum of the system. Thus, one finds 

To obtain the shear viscosity we have to know the value of h(p). For this purpose we linearize the transport equation 
Eq. (Plj) . The linearized collision term reads 



SC ^ 



2 



2EpT 

so that the linearized Boltzmann equation reads 



^ F,,[h{p)]V., , (40) 



J_,3i^_+AAp,^^V.^ - ^F^MpW., ■ (41) 
With this last result one realizes that the shear viscosity can also be written as 

^-^J ^P^Mp)Mh{p)] ■ (42) 

It is possible to solve the transport equation using variational methods. We will follow here the same strategy as 
that followed in Refs. [ill [ill, ^^^^ select the trial functions as 

oo 

Mp)=p"^6,i?,(p) , (43) 

where n is a parameter that will be determined in a variational procedure, and Bs{p) are orthogonal polynomials of 
order s and are defined such that the coefficient of the highest power is 1, and the orthogonality condition 

^^^^^^^^P^,P^,p-Br{p)BM = Ai-^Srs (44) 
is satisfied. 

Using this form of the solution one can check that Eq. (|39p gives 



where 



< = rTT+;rr(6 + «)C(5 + n) , (46) 



where T{z) and ^(2;) stand for the Gamma and Riemann zeta functions, respectively. 
On the other hand, the expression of the shear viscosity given by Eq. reads 



N=oo 

11= ^sbtMst , (47) 

s,t=0 
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where 



and 



Mst = J27:)^S^^Hp + K -P'- K')\\M\'f;Vl\l + + A^'^)A^A^ , (48) 



where we have defined the shorthand notation 



(50) 



Requiring that the two forms of the shear viscosity, Eqs. and (|T7|) . be equal implies 



^Ai-\M-'),o , (51) 



5T 



so that 



4ct 
25T2 



(4"V(M-^)oo . (52) 



In practical terms one performs the study by limiting the number of orthogonal polynomials included in the analysis. 
One can prove that 

^>^(4"^f(M-^)oo, (53) 

for a particular value of n and the number of orthogonal polynomials considered in the study. 

The evaluation of Eq. (|48|) requires the calculation of a five-dimensional integral. This can be realized after elim- 
inating from Eq. (|48p the p ' integral using the J-function of momentum conservation. The 5-function of energy 
conservation eliminates the integral over the magnitude of fc '. Then, we are left with a eight-dimensional integral 
that can be further simplified selecting the z axis along the vector p and realizing that only the difference in the 
two remaining azimuthal angles matters. Therefore, the problem is reduced to a five-dimensional integral over the 
magnitudes p and fc, two polar angles corresponding to k and fc', and one azimuthal angle (see Appendix C of Ref. [l^l , 
or [l^l ). This calculation can be performed numerically using the Vegas Monte Carlo algorithm [2l|. 

After a very simple dimensional analysis, it is very easy to check that the shear viscosity due to the binary collisions 
of superfluid phonons scales with the temperature as as it occurs in other superfluid systems. 

A variational numerical analysis can be performed once one specifies the choice of the EoS, and has the proper 
values of all the phonon self-couplings. However, independently of the specific choice of the EoS, one gets the same 
form of the solution for different superfluid systems |10l4l2| . One finds that the value of ry is maximized for the choice 
n = — 1, and for the lowest number of polynomials one can take in the expansion ((43|) . We will thus use this form of 
the solution to get the numerical values of the shear viscosity in /3-stable nuclear matter. 



B. Numerical values for the viscosity in superfluid neutron matter 

In this Subsection we present numerical results for the shear viscosity in /3-stable nuclear matter, using the APR 
EoS discussed in Sec. IIIII We have not been able to find a closed analytical expression for the value of 77, as the 
scattering cross section depends on all the phonon self- couplings g,rig, X,ri\,i and ?7a,2, and all diagrams in Fig. [T] 
contribute to the value of rj. 

We will present numerical results for values of the density higher and lower than 0.5 poj where po = 0.16 fm^'^ w 
2.8 10^^ g/(cm s) is the nuclear saturation density. It should be noted here that densities below 0.5 po correspond to 
the densities in the crust, while densities above that value would correspond to those of the core of the star. Our 
computation for the inner crust is valid for the superfluid phonon contribution to 77, but one should be aware that in 
the inner crust there are other light degrees of freedom that might have important contributions to the viscosity as 
well [l^j. We leave a complete study of 77 in the inner crust for future studies. 

In Fig.[3]we display a logarithmic plot of the shear viscosity where the expected 1/T^ dependence is clearly detected. 
There we present the results for four different values of the density: 0.25 po, 0.5 po, po and 2po. It is also clearly seen 
that 77 increases very rapidly when the density is increased. 



10 




Naively one could think that rj o: p^, traeing the direct dependence of the phonon self-couplings g and A on 
p. However, these couplings also depend on powers of the speed of sound and its derivatives, which are density 
dependent quantities, and whose values are different for different choices of the EoS. Therefore the p dependence of 
77 is not universal for all superfluids, contrary to what happens for the temperature dependence. 

By exploiting that rjT^ is a temperature-independent quantity, we study the p dependence of rj in Fig. 2] We 
consider both the APR EoS that describes /3-stable nuclear matter and the EoS of a free neutron gas. Amazingly, 
one observes that the viscosity increases by almost three orders of magnitude when going from po to 2 po for /3-stable 
matter, the increase being much more moderate for a free neutron gas. This fact is related to the behavior of Cs seen 
in Fig. [5] Our plot in Fig. 2] teaches us that the different choices of EoS have a very clear impact in the numerical 
values of ry. 

A relevant discussion also is to know the temperature regime where our results arc valid. The density of superfluid 
phonons becomes very dilute at very low T, and then it might be difficult to maintain a hydrodynamical description 
of their behavior. A relevant parameter that measures the validity of hydrodynamics is the mean free path. Hydro- 
dynamics is only valid when the mean free path is smaller that the typical macroscopic length of the system, in this 
case, the radius of the star. We show in Fig.[S]the mean free path of phonons within /3-stable matter for 0.5po, po and 
2po as a function of temperature. We also indicate an estimate of the limit of 10 Km, the radius of the star, above 
which phonon hydrodynamics should become invalid in neutron stars. The mean free path was extracted from the 
computation of 77 as [IJI 

I = — ^ , (54) 
n < p > 

where < p > is the thermal average momentum, and n is the phonon density: 



<p>^2.7-, „= / ^/, = ^(3)— ^. (55) 



The plots in Fig. [5] indicate that for temperatures below T ~ 0.1 MeV, the phonon mean free path is bigger than 
the size of the star. This would seem to indicate that for T < 0.1 MeV it is questionable to have a hydrodynamical 
description of the phonons of the star. We note that the critical temperature for the phase transition to the normal 
phase is Tc ~ 1 MeV, although this value, as well as the values we obtained for the phonon mean free path, are model 
dependent, that is, they depend on how one models the nucleonic interactions and the EoS of the system. We should 
also warn the reader that we have only made a tentative estimate, as the phonons also collide with other particles in 
the star, such as the electrons. To assess the mean free path of the phonons due to their collision with the electrons, 
one should derive an effective field theory describing their interactions. We leave this computation for future studies. 
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FIG. 4: Shear viscosity multiplied by irj T^) as function of density for a free neutron gas and /3-stabIe nuclear matter. We 
also plot rj due to electron-electron scattering, according to the estimate of Cutler and Lindblom, see Ref [i]. As the last 
is a T-dependent quantity, we chose T = 0.1 MeV in that curve. 



P-stable matter 



le+16 




T[MeV] 

FIG. 5: Phonon mean-free path in /3-stable matter as function of temperature for different densities. The straight line indicates 
the hydrodynamic limit: Ih < 10 Km. 



C. Comparison with the shear viscosity due to phonons of other superfluid systems 



Having set a general formulation for the computation of the shear viscosity due to the binary collisions of superfluid 
phonons, we would like to compare here our results with those obtained in the literature for other superfluid systems. 

Landau and Khalatnikov were the first to compute the viscosity due to the phonons in the regime where they 
were supposed to dominate the transport phenomena in the cold regime of *He [2^. Landau postulated an ad hoc 
Hamiltonian, valid for ^He. The formalism we have followed to model the phonon interactions is very different than 
the one used for ^He, but as shown in Ref. Q, after a Legendre transformation one can convert the Lagrangian Eq. ([T]) 
into the Hamiltonian used by Landau and Khalatnikov in their computations. Thus, one would expect to have very 
similar results for the viscosity. For T < 0.8 K, when the phonons dominate. Landau and Khalatnikov obtained \2c 



4^ 



2 2 

P c,- 



5 • 2 



17 



3-13! 



1 j_ _P_§£s. 



rp5 



(56) 



However, this result was obtained after making several approximations. More particularly, after considering that the 
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fourth order term in a perturbative expansion of the Hamiltonian is ncghgiblc (in our language, the contact interaction 
diagram), and considering only collisions where the momentum of one of the particle that collides is much smaller 
than the momentum of the other one. While the scaling behavior with the temperature T is the same as the one we 
found out, the numerical prefactor is not. The reason is due to the approximations that Landau and Khalatnikov 
performed to obtain the scattering cross section. In fact, our numerical analysis suggests that in our superfluid all 
the different Feynman diagrams contribute to 77, and that also not only collisions where p <^ k are relevant in the 
computation of 77. 

With our computations wc can easily get the value of the shear viscosity for the ultracold Fermi gas in the unitarity 
limit [111 . The EoS for that system is Q 

_ 4V2m3^ 5/2 .57^ 

where ^ is the so-called Bertsch parameter, which relates the chemical potential and the Fermi energy fi = ^Ep- One 
then easily finds the value of p and also that 

^2^2^ P_dc^^l fPd^^_2 ^^^^ 

^ 3m ' Cs dp 3 ' Cs dp^ 9 

The computation of the shear viscosity for the superfluid phase was performed by Rupak and Schafcr. They obtained 

^ = 9.3. 10-6 -| , (59) 
where Tp ^ fi/^. We reproduce the same scaling behavior, but find a small discrepancy in the numerical pre- factor: 

= 3-4 • 10-4^- (60) 



c 



Because the computations in Ref. [llj and ours are essentially the same, the numerical disagreement seems to be 
related to a different pre-factor of the scattering cross section and the imaginary part of the phonon self-energy. 



V. DISCUSSION 



We have set a formulation of the binary collisions of the phonons and compute the associated shear viscosity of 
a non-relativistic superfluid in terms of the equation of state of the system. Our computations are rather general, 
and can be used for different superfluid systems, provided they share the same underlying symmetries f26| . We have 
exploited the strength and universal character of the effective field theory that describes the phonon interactions 
at leading order in a momentum and energy expansion. In this way, we see that the binary collisions of phonons 
produce a shear viscosity that scales with the temperature as 7] cx l/T^, and this is universal for all superfluids sharing 
the same symmetries. The density dependence of i] is not universal, and depends on the EoS of the system under 
consideration. 

Our primary goal in this manuscript was computing the phonon contribution to the shear viscosity in superfluid 
neutron stars, if the neutrons pair in a channel. For the superfluid in the inner crust and for neutron pairing in 
a channel, other contributions to tj are as well important and will be studied somewhere else. 

We believe that the results we obtained here should be considered in the study of the behavior of different oscillation 
modes of superfluid neutron stars. While in the normal phase the viscosity is dominated by neutron-neutron scattering, 
in the superfluid phase it was considered that the viscosity should be dominated by electron-electron scattering. There 
are several estimates of the electron contribution to the shear viscosity (see Fig. 1 of Ref. Q for a comparison). We 
note that the phonon contribution to rj is several orders of magnitude higher than the electron contribution in a 
certain low temperature range (see Fig. |4] for an explicit comparison of T^r] at T = 0.1 MeV) , but this statement also 
depends on how deep inside one looks in the star, the phonon viscosity being much larger when the density increases. 
It should be also important to assess when phonons are in a ballistic regime or not. In any case, our results indicate 
that the effect of phonons in the value of the transport coefficients cannot be ignored, and that one should check their 
impact in the study of the oscillation modes of a superfluid neutron star. 
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